Clinging to the top: natal dispersal tracks climate gradient in a trailing-edge population of a migratory songbird.

PURPOSE: Trailing-edge populations at the low-latitude, receding edge of a shifting range face high extinction risk from climate change unless they are able to track optimal environmental conditions through dispersal. METHODS: We fit dispersal models to the locations of 3165 individually-marked black-throated blue warblers (Setophaga caerulescens) in the southern Appalachian Mountains in North Carolina, USA from 2002 to 2023. Black-throated blue warbler breeding abundance in this population has remained relatively stable at colder and wetter areas at higher elevations but has declined at warmer and drier areas at lower elevations. RESULTS: Median dispersal distance of young warblers was 917 m (range 23-3200 m), and dispersal tended to be directed away from warm and dry locations. In contrast, adults exhibited strong site fidelity between breeding seasons and rarely dispersed more than 100 m (range 10-1300 m). Consequently, adult dispersal kernels were much more compact and symmetric than natal dispersal kernels, suggesting adult dispersal is unlikely a driving force of declines in this population. CONCLUSION: Our findings suggest that directional natal dispersal may mitigate fitness costs for trailing-edge populations by allowing individuals to track changing climate and avoid warming conditions at warm-edge range boundaries.

Ecological filtering shapes the impacts of agricultural deforestation on biodiversity.

The biodiversity impacts of agricultural deforestation vary widely across regions. Previous efforts to explain this variation have focused exclusively on the landscape features and management regimes of agricultural systems, neglecting the potentially critical role of ecological filtering in shaping deforestation tolerance of extant species assemblages at large geographical scales via selection for functional traits. Here we provide a large-scale test of this role using a global database of species abundance ratios between matched agricultural and native forest sites that comprises 71 avian assemblages reported in 44 primary studies, and a companion database of 10 functional traits for all 2,647 species involved. Using meta-analytic, phylogenetic and multivariate methods, we show that beyond agricultural features, filtering by the extent of natural environmental variability and the severity of historical anthropogenic deforestation shapes the varying deforestation impacts across species assemblages. For assemblages under greater environmental variability-proxied by drier and more seasonal climates under a greater disturbance regime-and longer deforestation histories, filtering has attenuated the negative impacts of current deforestation by selecting for functional traits linked to stronger deforestation tolerance. Our study provides a previously largely missing piece of knowledge in understanding and managing the biodiversity consequences of deforestation by agricultural deforestation.

Genomics-informed conservation units reveal spatial variation in climate vulnerability in a migratory bird.

Identifying genetic conservation units (CUs) in threatened species is critical for the preservation of adaptive capacity and evolutionary potential in the face of climate change. However, delineating CUs in highly mobile species remains a challenge due to high rates of gene flow and genetic signatures of isolation by distance. Even when CUs are delineated in highly mobile species, the CUs often lack key biological information about what populations have the most conservation need to guide management decisions. Here we implement a framework for CU identification in the Canada Warbler (Cardellina canadensis), a migratory bird species of conservation concern, and then integrate demographic modelling and genomic offset to guide conservation decisions. We find that patterns of whole genome genetic variation in this highly mobile species are primarily driven by putative adaptive variation. Identification of CUs across the breeding range revealed that Canada Warblers fall into two evolutionarily significant units (ESU), and three putative adaptive units (AUs) in the South, East, and Northwest. Quantification of genomic offset, a metric of genetic changes necessary to maintain current gene-environment relationships, revealed significant spatial variation in climate vulnerability, with the Northwestern AU being identified as the most vulnerable to future climate change. Alternatively, quantification of past population trends within each AU revealed the steepest population declines have occurred within the Eastern AU. Overall, we illustrate that genomics-informed CUs provide a strong foundation for identifying current and future regional threats that can be used to inform management strategies for a highly mobile species in a rapidly changing world.

Climate-mediated population dynamics of a migratory songbird differ between the trailing edge and range core.

Understanding the demographic drivers of range contractions is important for predicting species' responses to climate change; however, few studies have examined the effects of climate change on survival and recruitment across species' ranges. We show that climate change can drive trailing edge range contractions through the effects on apparent survival, and potentially recruitment, in a migratory songbird. We assessed the demographic drivers of trailing edge range contractions using a long-term demography dataset for the black-throated blue warbler (Setophaga caerulescens) collected across elevational climate gradients at the trailing edge and core of the breeding range. We used a Bayesian hierarchical model to estimate the effect of climate change on apparent survival and recruitment and to forecast population viability at study plots through 2040. The trailing edge population at the low-elevation plot became locally extinct by 2017. The local population at the mid-elevation plot at the trailing edge gradually declined and is predicted to become extirpated by 2040. Population declines were associated with warming temperatures at the mid-elevation plot, although results were more equivocal at the low-elevation plot where we had fewer years of data. Population density was stable or increasing at the range core, although warming temperatures are predicted to cause population declines by 2040 at the low-elevation plot. This result suggests that even populations within the geographic core of the range are vulnerable to climate change. The demographic drivers of local population declines varied between study plots, but warming temperatures were frequently associated with declining rates of population growth and apparent survival. Declining apparent survival in our study system is likely to be associated with increased adult emigration away from poor-quality habitats. Our results suggest that demographic responses to warming temperatures are complex and dependent on local conditions and geographic range position, but spatial variation in population declines is consistent with the climate-mediated range shift hypothesis. Local populations of black-throated blue warblers near the warm-edge range boundary at low latitudes and low elevations are likely to be the most vulnerable to climate change, potentially leading to local extirpation and range contractions.

Monitoring partially marked populations using camera and telemetry data.

Long-term monitoring is an important component of effective wildlife conservation. However, many methods for estimating density are too costly or difficult to implement over large spatial and temporal extents. Recently developed spatial mark-resight (SMR) models are increasingly being applied as a cost-effective method to estimate density when data include detections of both marked and unmarked individuals. We developed a generalized SMR model that can accommodate long-term camera data and auxiliary telemetry data for improved spatiotemporal inference in monitoring efforts. The model can be applied in two stages, with detection parameters estimated in the first stage using telemetry data and camera detections of instrumented individuals. Density is estimated in the second stage using camera data, with all individuals treated as unmarked. Serial correlation in detection and density parameters is accounted for using time-series models. The two-stage approach reduces computational demands and facilitates the application to large data sets from long-term monitoring initiatives. We applied the model to 3 years (2015-2017) of white-tailed deer (Odocoileus virginianus) data collected in three study areas of the Big Cypress Basin, Florida, USA. In total, 59 females marked with ear tags and fitted with GPS-telemetry collars were detected along with unmarked females on 180 remote cameras. Most of the temporal variation in density was driven by seasonal fluctuations, but one study area exhibited a slight population decline during the monitoring period. Modern technologies such as camera traps provide novel possibilities for long-term monitoring, but the resulting massive data sets, which are subject to unique sources of observation error, have posed analytical challenges. The two-stage spatial mark-resight framework provides a solution with lower computational demands than joint SMR models, allowing for easier implementation in practice. In addition, after detection parameters have been estimated, the model may be used to estimate density even if no synchronous auxiliary information on marked individuals is available, which is often the case in long-term monitoring.

An integrated path for spatial capture-recapture and animal movement modeling.

Ecologists and conservation biologists increasingly rely on spatial capture-recapture (SCR) and movement modeling to study animal populations. Historically, SCR has focused on population-level processes (e.g., vital rates, abundance, density, and distribution), whereas animal movement modeling has focused on the behavior of individuals (e.g., activity budgets, resource selection, migration). Even though animal movement is clearly a driver of population-level patterns and dynamics, technical and conceptual developments to date have not forged a firm link between the two fields. Instead, movement modeling has typically focused on the individual level without providing a coherent scaling from individual- to population-level processes, whereas SCR has typically focused on the population level while greatly simplifying the movement processes that give rise to the observations underlying these models. In our view, the integration of SCR and animal movement modeling has tremendous potential for allowing ecologists to scale up from individuals to populations and advancing the types of inferences that can be made at the intersection of population, movement, and landscape ecology. Properly accounting for complex animal movement processes can also potentially reduce bias in estimators of population-level parameters, thereby improving inferences that are critical for species conservation and management. This introductory article to the Special Feature reviews recent advances in SCR and animal movement modeling, establishes a common notation, highlights potential advantages of linking individual-level (Lagrangian) movements to population-level (Eulerian) processes, and outlines a general conceptual framework for the integration of movement and SCR models. We then identify important avenues for future research, including key challenges and potential pitfalls in the developments and applications that lie ahead.

Modeling abundance, distribution, movement and space use with camera and telemetry data.

Studies of animal abundance and distribution are often conducted independently of research on movement, despite the important links between processes. Movement can cause rapid changes in spatial variation in density, and movement influences detection probability and therefore estimates of abundance from inferential methods such as spatial capture-recapture (SCR). Technological developments including camera traps and GPS telemetry have opened new opportunities for studying animal demography and movement, yet statistical models for these two data types have largely developed along parallel tracks. We present a hierarchical model in which both datasets are conditioned on a movement process for a clearly defined population. We fitted the model to data from 60 camera traps and 23,572 GPS telemetry locations collected on 17 male white-tailed deer in the Big Cypress National Preserve, Florida, USA during July 2015. Telemetry data were collected on a 3-4 h acquisition schedule, and we modeled the movement paths of all individuals in the region with a Ornstein-Uhlenbeck process that included individual-specific random effects. Two of the 17 deer with GPS collars were detected on cameras. An additional 20 male deer without collars were detected on cameras and individually identified based on their unique antler characteristics. Abundance was 126 (95% CI: 88-177) in the 228 km2 region, only slightly higher than estimated using a standard SCR model: 119 (84-168). The standard SCR model, however, was unable to describe individual heterogeneity in movement rates and space use as revealed by the joint model. Joint modeling allowed the telemetry data to inform the movement model and the SCR encounter model, while leveraging information in the camera data to inform abundance, distribution and movement. Unlike most existing methods for population-level inference on movement, the joint SCR-movement model can yield unbiased inferences even if non-uniform sampling is used to deploy transmitters. Potential extensions of the model include the addition of resource selection parameters, and relaxation of the closure assumption when interest lies in survival and recruitment. These developments would contribute to the emerging holistic framework for the study of animal ecology, one that uses modern technology and spatio-temporal statistics to learn about interactions between behavior and demography.

Effects of Bait on Male White-Tailed Deer Resource Selection.

Bait is often used to increase wildlife harvest susceptibility, enhance viewing opportunities, and survey wildlife populations. The effects of baiting depend on how bait influences space use and resource selection at multiple spatial scales. Although telemetry studies allow for inferences about resource selection within home ranges (third-order selection), they provide limited information about spatial variation in density, which is the result of second-order selection. Recent advances in spatial capture-recapture (SCR) techniques allow exploration of second- and third-order selection simultaneously using non-invasive methods such as camera traps. Our objectives were to describe how short-term baiting affects white-tailed deer (Odocoileus virginianus) behavior and distribution. We fit SCR models to camera data from baited and unbaited locations in southwestern Georgia to assess the effects of short-term baiting on second- and third-order selection of deer during summer and winter surveys. We found little evidence of second-order selection during late summer or early winter surveys when camera surveys using bait are typically conducted. However, we found evidence for third-order selection, indicating that resource selection within home ranges is affected. Concentrations in space use resulting from baiting may enhance disease transmission, change harvest susceptibility, and potentially bias the outcome of camera surveys using bait.

Optimizing conservation in species-specific agricultural landscapes.

Recovery of grassland birds in agricultural landscapes is a global imperative. Agricultural landscapes are complex, and the value of resource patches may vary substantially among species. The spatial extent at which landscape features affect populations (i.e., scale of effect) may also differ among species. There is a need for regional-scale conservation planning that considers landscape-scale and species-specific responses of grassland birds to environmental change. We developed a spatially explicit approach to optimizing grassland conservation in the context of species-specific landscapes and prioritization of species recovery and applied it to a conservation program in Kentucky (USA). We used a hierarchical distance-sampling model with an embedded scale of effect predictor to estimate the relationship between landscape structure and abundance of eastern meadowlarks (Sturnella magna), field sparrows (Spizella pusilla), and northern bobwhites (Colinus virginianus). We used a novel spatially explicit optimization procedure rooted in multi-attribute utility theory to design alternative conservation strategies (e.g., prioritize only northern bobwhite recovery or assign equal weight to each species' recovery). Eastern meadowlarks and field sparrows were more likely to respond to landscape-scale resource patch adjacencies than landscape-scale patch densities. Northern bobwhite responded to both landscape-scale resource patch adjacencies and densities and responded strongly to increased grassland density. Effects of landscape features on local abundance decreased as distance increased and had negligible influence at 0.8 km for eastern meadowlarks (0.7-1.2 km 95% Bayesian credibility intervals [BCI]), 2.5 km for field sparrows (1.5-5.8 km 95% BCI), and 8.4 km for bobwhite (6.4-26 km 95% BCI). Northern bobwhites were predicted to benefit greatly from future grassland conservation regardless of conservation priorities, but eastern meadowlark and field sparrow were not. Our results suggest similar species can respond differently to broad-scale conservation practices because of species-specific, distance-dependent relationships with landscape structure. Our framework is quantitative, conceptually simple, customizable, and predictive and can be used to optimize conservation in heterogeneous ecosystems while considering landscape-scale processes and explicit prioritization of species recovery.

La recuperación de las aves de pastizal en los paisajes agrícolas es una obligación mundial. Los paisajes agrícolas son complejos y el valor de los fragmentos con recursos puede variar sustancialmente entre especies. La magnitud espacial a la que las características del paisaje afectan a las poblaciones (es decir, la escala del efecto) también puede diferir entre especies. Existe la necesidad de una planeación de la conservación a escala regional que considere la escala de paisaje y las respuestas específicas de especie de aves de pastizal al cambio ambiental. Desarrollamos una estrategia espacialmente explícita para optimizar la conservación de pastizales en el contexto de los paisajes de especies específicas y la priorización de la recuperación de especies y la aplicamos a un programa de conservación en Kentucky (E.U.A.). Usamos un modelo jerárquico de muestreo a distancia con una escala integrada del efecto pronosticador para estimar la relación entre la estructura del paisaje y la abundancia de la alondra oriental de pradera (Sturnella magna), el gorrión de campo (Spizella pusilla) y la codorniz norteña (Colinus virginianus). Usamos un novedoso procedimiento de optimización espacialmente explícito basado en la teoría de utilidad multicaracterística para diseñar estrategias de conservación alternativas (p. ej.: priorizar solamente la recuperación de la codorniz norteña o asignar una importancia idéntica a la recuperación de cada especie). La alondra y el gorrión tuvieron una mayor probabilidad de responder a la proximidad de fragmentos con recursos a escala de paisaje que a la densidad de fragmentos a escala de paisaje. La codorniz respondió tanto a la proximidad de fragmentos con recursos a escala de paisaje como a la densidad y también respondió fuertemente al incremento en la densidad del pastizal. Los efectos de las características del paisaje sobre la abundancia local disminuyeron conforme incrementó la distancia, representando una influencia insignificante a los 0.8 km para la alondra (0.7-1.2 km 95% de intervalos de credibilidad bayesiana [ICB]), a los 2.5 km para el gorrión (1.5-5.8 km 95% ICB) y a los 8.4 km para la codorniz (6.4-26 km 95% ICB). Se pronosticó que la codorniz se beneficiaría enormemente con la conservación futura de los pastizales sin importar las prioridades de conservación, pero no fue el caso para la alondra y el gorrión. Nuestros resultados sugieren que especies similares pueden responder de manera diferente a las prácticas de conservación a escalas generalizadas debido a las relaciones específicas de especie y dependientes de la distancia con la estructura del paisaje. Nuestro marco de trabajo es cuantitativo, conceptualmente simple, adaptable y predictivo y puede usarse para optimizar la conservación en los ecosistemas heterogéneos a la vez que considera los procesos a escala de paisaje y la priorización explícita de la recuperación de las especies.

An experimental test of the Allee effect range limitation hypothesis.

Understanding how climate change impacts trailing-edge populations requires information about how abiotic and biotic factors limit their distributions. Theory indicates that socially mediated Allee effects can limit species distributions by suppressing growth rates of peripheral populations when social information is scarce. The goal of our research was to determine if socially mediated Allee effects limit the distribution of Canada warbler Cardellina canadensis at the trailing-edge of the geographic range. Using 4 years of observational data from 71 sites and experimental data at 10 sites, we tested two predictions of the socially mediated range limitation hypothesis: (a) local growth rates should be positively correlated with local density and (b) the addition of social cues immediately outside the trailing-edge range boundary would result in colonization of formerly unoccupied habitat and increased growth rates. During the third breeding season, social cues were experimentally added at 10 formerly unoccupied sites within and beyond the species' local range margin to determine if the addition of social information could increase density and effectively expand the species' range. No experimental sites were colonized after adding social cues and no evidence of Allee effects was found. Rather, temperature, precipitation and negative density dependence strongly influenced population growth rates. Although theoretical models indicate that the presence of socially mediated Allee effects at species range boundaries could increase the rate of climate-induced range shifts and local extinctions, empirical results from the first test of this hypothesis suggest that Allee effects play a minimal role in limiting species' distributions.

A statistical forecasting approach to metapopulation viability analysis.

Conservation of at-risk species is aided by reliable forecasts of the consequences of environmental change and management actions on population viability. Forecasts from conventional population viability analysis (PVA) are made using a two-step procedure in which parameters are estimated, or elicited from expert opinion, and then plugged into a stochastic population model without accounting for parameter uncertainty. Recently developed statistical PVAs differ because forecasts are made conditional on models fitted to empirical data. The statistical forecasting approach allows for uncertainty about parameters, but it has rarely been applied in metapopulation contexts where spatially explicit inference is needed about colonization and extinction dynamics and other forms of stochasticity that influence metapopulation viability. We conducted a statistical metapopulation viability analysis (MPVA) using 11 yr of data on the federally threatened Chiricahua leopard frog (Lithobates chiricahuensis) to forecast responses to landscape heterogeneity, drought, environmental stochasticity, and management. We evaluated several future environmental scenarios and pond restoration options designed to reduce extinction risk. Forecasts over a 50-yr time horizon indicated that metapopulation extinction risk was <4% for all scenarios, but uncertainty was high. Without pond restoration, extinction risk is forecasted to be 3.9% (95% CI 0-37%) by year 2066. Restoring six ponds by increasing their hydroperiod reduced extinction risk to <1% and greatly reduced uncertainty (95% CI 0-2%). Our results suggest that managers can mitigate the impacts of drought and environmental stochasticity on metapopulation viability by maintaining ponds that hold water throughout the year and keeping them free of invasive predators. Our study illustrates the utility of the spatially explicit statistical forecasting approach to MPVA in conservation planning efforts.

Chronology of reproductive investment determines predation risk aversion in a felid-ungulate system.

Fear of predators can behaviorally mediate prey population dynamics, particularly when predation risk influences reproductive investment. However, the costs of reproductive investment may mitigate predation risk aversion relative to periods when the link between reproductive output and prey behavior is weaker.We posit that intensity of reproductive investment in ungulates may predict their response to predation risk such that the sexes increase risk exposure during biological seasons that are pivotal to reproductive success, such as the fawn-rearing and breeding seasons for females and males, respectively.We examined the activity patterns of sympatric white-tailed deer (Odocoileus virginianus), a sexually segregated polygynous ungulate, and Florida panthers (Puma concolor coryi) in the context of the "risky times - risky places hypothesis" and the reproductive strategy hypothesis. We compared detection rates and diel activity overlap of both species using motion-triggered camera traps positioned on (n = 120) and off (n = 60) anthropogenic trails across five reproductive seasons.Florida panthers were nocturnal and primarily observed on-trail providing an experimental framework with risky times and risky places. Contrary to studies in other taxa inversely correlating prey reproductive investment to predation risk, the sexes of deer were more risk prone during sex-specific seasons associated with intense reproductive investment.Our results suggest spatiotemporally variable predation risk influences sex-specific behavioral decision-making in deer such that reproductive success is maximized.

Altitudinal migration: ecological drivers, knowledge gaps, and conservation implications.

Animal migration has been the subject of intensive research for more than a century, but most research has focused on long-distance rather than short-distance migration. Altitudinal migration is a form of short-distance migration in which individuals perform seasonal elevational movements. Despite its geographic and taxonomic ubiquity, there is relatively little information about the intrinsic and extrinsic factors that influence altitudinal migratory behaviour. Without this information, it is difficult to predict how rapid environmental changes will affect population viability of altitudinal migrants. To synthesize current knowledge, we compiled literature on altitudinal migration for all studied taxa, and identified the leading hypotheses explaining this behaviour. Studies of animal altitudinal migration cover many taxonomic lineages, with birds being the most commonly studied group. Altitudinal migration occurs in all continents except for Antarctica, but about a third of the literature focused on altitudinal migration in North America. Most research suggests that food and weather are the primary extrinsic drivers of altitudinal migration. In addition, substantial individual-level variation in migratory propensity exists. Individual characteristics that are associated with sex, dominance rank, and body size explain much of the variation in migratory propensity in partially migratory populations, but individual-level correlates are poorly known for most taxa. More research is needed to quantify the effects of habitat loss, habitat fragmentation, and climate change on altitudinal migrants. Demographic studies of individually marked populations would be particularly valuable for advancing knowledge of the cascading effects of environmental change on migratory propensity, movement patterns, and population viability. We conclude our review with recommendations for study designs and modelling approaches that could be used to narrow existing knowledge gaps, which currently hinder effective conservation of altitudinal migratory species.

Increasing connectivity between metapopulation ecology and landscape ecology.

Metapopulation ecology and landscape ecology aim to understand how spatial structure influences ecological processes, yet these disciplines address the problem using fundamentally different modeling approaches. Metapopulation models describe how the spatial distribution of patches affects colonization and extinction, but often do not account for the heterogeneity in the landscape between patches. Models in landscape ecology use detailed descriptions of landscape structure, but often without considering colonization and extinction dynamics. We present a novel spatially explicit modeling framework for narrowing the divide between these disciplines to advance understanding of the effects of landscape structure on metapopulation dynamics. Unlike previous efforts, this framework allows for statistical inference on landscape resistance to colonization using empirical data. We demonstrate the approach using 11 yr of data on a threatened amphibian in a desert ecosystem. Occupancy data for Lithobates chiricahuensis (Chiricahua leopard frog) were collected on the Buenos Aires National Wildlife Refuge (BANWR), Arizona, USA from 2007 to 2017 following a reintroduction in 2003. Results indicated that colonization dynamics were influenced by both patch characteristics and landscape structure. Landscape resistance increased with increasing elevation and distance to the nearest streambed. Colonization rate was also influenced by patch quality, with semi-permanent and permanent ponds contributing substantially more to the colonization of neighboring ponds relative to intermittent ponds. Ponds that only hold water intermittently also had the highest extinction rate. Our modeling framework can be widely applied to understand metapopulation dynamics in complex landscapes, particularly in systems in which the environment between habitat patches influences the colonization process.

Estimating population extinction thresholds with categorical classification trees for Louisiana black bears.

Monitoring vulnerable species is critical for their conservation. Thresholds or tipping points are commonly used to indicate when populations become vulnerable to extinction and to trigger changes in conservation actions. However, quantitative methods to determine such thresholds have not been well explored. The Louisiana black bear (Ursus americanus luteolus) was removed from the list of threatened and endangered species under the U.S. Endangered Species Act in 2016 and our objectives were to determine the most appropriate parameters and thresholds for monitoring and management action. Capture mark recapture (CMR) data from 2006 to 2012 were used to estimate population parameters and variances. We used stochastic population simulations and conditional classification trees to identify demographic rates for monitoring that would be most indicative of heighted extinction risk. We then identified thresholds that would be reliable predictors of population viability. Conditional classification trees indicated that annual apparent survival rates for adult females averaged over 5 years ([Formula: see text]) was the best predictor of population persistence. Specifically, population persistence was estimated to be ≥95% over 100 years when [Formula: see text], suggesting that this statistic can be used as threshold to trigger management intervention. Our evaluation produced monitoring protocols that reliably predicted population persistence and was cost-effective. We conclude that population projections and conditional classification trees can be valuable tools for identifying extinction thresholds used in monitoring programs.

An open-population hierarchical distance sampling model.

Modeling population dynamics while accounting for imperfect detection is essential to monitoring programs. Distance sampling allows estimating population size while accounting for imperfect detection, but existing methods do not allow for estimation of demographic parameters. We develop a model that uses temporal correlation in abundance arising from underlying population dynamics to estimate demographic parameters from repeated distance sampling surveys. Using a simulation study motivated by designing a monitoring program for Island Scrub-Jays (Aphelocoma insularis), we investigated the power of this model to detect population trends. We generated temporally autocorrelated abundance and distance sampling data over six surveys, using population rates of change of 0.95 and 0.90. We fit the data generating Markovian model and a mis-specified model with a log-linear time effect on abundance, and derived post hoc trend estimates from a model estimating abundance for each survey separately. We performed these analyses for varying numbers of survey points. Power to detect population changes was consistently greater under the Markov model than under the alternatives, particularly for reduced numbers of survey points. The model can readily be extended to more complex demographic processes than considered in our simulations. This novel framework can be widely adopted for wildlife population monitoring.

Modeling structured population dynamics using data from unmarked individuals.

The study of population dynamics requires unbiased, precise estimates of abundance and vital rates that account for the demographic structure inherent in all wildlife and plant populations. Traditionally, these estimates have only been available through approaches that rely on intensive mark-recapture data. We extended recently developed N-mixture models to demonstrate how demographic parameters and abundance can be estimated for structured populations using only stage-structured count data. Our modeling framework can be used to make reliable inferences on abundance as well as recruitment, immigration, stage-specific survival, and detection rates during sampling. We present a range of simulations to illustrate the data requirements, including the number of years and locations necessary for accurate and precise parameter estimates. We apply our modeling framework to a population of northern dusky salamanders (Desmognathus fuscus) in the mid-Atlantic region (USA) and find that the population is unexpectedly declining. Our approach represents a valuable advance in the estimation of population dynamics using multistate data from unmarked individuals and should additionally be useful in the development of integrated models that combine data from intensive (e.g., mark-recapture) and extensive (e.g., counts) data sources.

Inferences about population dynamics from count data using multistate models: a comparison to capture-recapture approaches.

Wildlife populations consist of individuals that contribute disproportionately to growth and viability. Understanding a population's spatial and temporal dynamics requires estimates of abundance and demographic rates that account for this heterogeneity. Estimating these quantities can be difficult, requiring years of intensive data collection. Often, this is accomplished through the capture and recapture of individual animals, which is generally only feasible at a limited number of locations. In contrast, N-mixture models allow for the estimation of abundance, and spatial variation in abundance, from count data alone. We extend recently developed multistate, open population N-mixture models, which can additionally estimate demographic rates based on an organism's life history characteristics. In our extension, we develop an approach to account for the case where not all individuals can be assigned to a state during sampling. Using only state-specific count data, we show how our model can be used to estimate local population abundance, as well as density-dependent recruitment rates and state-specific survival. We apply our model to a population of black-throated blue warblers (Setophaga caerulescens) that have been surveyed for 25 years on their breeding grounds at the Hubbard Brook Experimental Forest in New Hampshire, USA. The intensive data collection efforts allow us to compare our estimates to estimates derived from capture-recapture data. Our model performed well in estimating population abundance and density-dependent rates of annual recruitment/immigration. Estimates of local carrying capacity and per capita recruitment of yearlings were consistent with those published in other studies. However, our model moderately underestimated annual survival probability of yearling and adult females and severely underestimates survival probabilities for both of these male stages. The most accurate and precise estimates will necessarily require some amount of intensive data collection efforts (such as capture-recapture). Integrated population models that combine data from both intensive and extensive sources are likely to be the most efficient approach for estimating demographic rates at large spatial and temporal scales.

Spatial capture--recapture models for jointly estimating population density and landscape connectivity.

Population size and landscape connectivity are key determinants of population viability, yet no methods exist for simultaneously estimating density and connectivity parameters. Recently developed spatial capture--recapture (SCR) models provide a framework for estimating density of animal populations but thus far have not been used to study connectivity. Rather, all applications of SCR models have used encounter probability models based on the Euclidean distance between traps and animal activity centers, which implies that home ranges are stationary, symmetric, and unaffected by landscape structure. In this paper we devise encounter probability models based on "ecological distance," i.e., the least-cost path between traps and activity centers, which is a function of both Euclidean distance and animal movement behavior in resistant landscapes. We integrate least-cost path models into a likelihood-based estimation scheme for spatial capture-recapture models in order to estimate population density and parameters of the least-cost encounter probability model. Therefore, it is possible to make explicit inferences about animal density, distribution, and landscape connectivity as it relates to animal movement from standard capture-recapture data. Furthermore, a simulation study demonstrated that ignoring landscape connectivity can result in negatively biased density estimators under the naive SCR model.

A small-scale land-sparing approach to conserving biological diversity in tropical agricultural landscapes.

Two contrasting strategies have been proposed for conserving biological diversity while meeting the increasing demand for agricultural products: land sparing and land sharing production systems. Land sparing involves increasing yield to reduce the amount of land needed for agriculture, whereas land-sharing agricultural practices incorporate elements of native ecosystems into the production system itself. Although the conservation value of these systems has been extensively debated, empirical studies are lacking. We compared bird communities in shade coffee, a widely practiced land-sharing system in which shade trees are maintained within the coffee plantation, with bird communities in a novel, small-scale, land-sparing coffee-production system (integrated open canopy or IOC coffee) in which farmers obtain higher yields under little or no shade while conserving an area of forest equal to the area under cultivation. Species richness and diversity of forest-dependent birds were higher in the IOC coffee farms than in the shade coffee farms, and community composition was more similar between IOC coffee and primary forest than between shade coffee and primary forest. Our study represents the first empirical comparison of well-defined land sparing and land sharing production systems. Because IOC coffee farms can be established by allowing forest to regenerate on degraded land, widespread adoption of this system could lead to substantial increases in forest cover and carbon sequestration without compromising agricultural yield or threatening the livelihoods of traditional small farmers. However, we studied small farms (<5 ha); thus, our results may not generalize to large-scale land-sharing systems. Furthermore, rather than concluding that land sparing is generally superior to land sharing, we suggest that the optimal approach depends on the crop, local climate, and existing land-use patterns.